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ABSTRACT 

We present a modern ephemeris and modern light curve of the first-discovered 
py/| Cepheid variable in M31, Edwin Hubble's M31-V1. Observers of the American 

Association of Variable Star Observers undertook these observations during the 
i-C ■ latter half of 2010. The observations were in support of an outreach program by 

the Space Telescope Science Institute's Hubble Heritage project, but the resulting 
2 ■ data are the first concentrated observations of M31-V1 made in modern times. 

AAVSO observers obtained 214 V-band, Rc-band, and unfiltered observations 
from which a current ephemeris was derived. The ephemeris derived from these 
observations is JD Max = 2455430.5(±0.5) + 31.4(±0.1)£. The period derived 
^ . from the 2010 data are in agreement with the historic values of the period, but 

the single season of data precludes a more precise determination of the period 
^ . or measurement of the period change using these data alone. However, using 

an ephemeris based upon the period derived by Baade and Swope we are able 
to fit all of the observed data acceptably well. Continued observations in the 
modern era will be very valuable in linking these modern data with data from 
the 1920s-30s and 1950s, and will enable us to measure period change in this 
historic Cepheid. In particular, we strongly encourage intensive observations of 
this star around predicted times of maximum to constrain the date of maximum 
to better than 0.5 days. 



Subject headings: Stars: variables: Cepheids - Stars: individual: M31-V1 - 
Galaxies: Local Group - Galaxies: individual: M31 - Astronomical databases: 
AAVSO - Facilities: AAVSO 



1 also Robert Ferguson Observatory, Sugarloaf Ridge State Park, Kenwood, California 
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Introduction 



The star M31-V1 (RA: 00:41: 27.30, Dec: +41:10 :10.4, J2000) was the first Cepheid dis- 
covered by Edwin Hubble in M31 (jChristiansonlll995l ). thus making i t an o bject of historical 
as well as astrophysical interest. The object was listed in Hubble (jl925 ) with a period of 
31 .41 days and a maximum magnitude of 18.2. Hubble published his light curves in full 



m 



Hubble! (119291 ). clearly demonstrating that this star and several others like it in M31 ex- 
hibited Cepheid-like behavior. After this point, M31 was established as being a unique and 
distant Galaxy, independent of the Milky Way. Little or no work was done o n these Cepheids 
until the work of W. Baade, published posthumously by H. Swope in 1965 ( jBaade fc Swope 
19651 ). Baade provided several additional years of data for M31-V1 along with many other 
variables. Baade's observations confirmed those of Hubble, and their analysis showed that 
the gross properties of M31-V1 were not significantly changed since the 1920s. Beyond these 
two data sets, little or no published time series photometry exists for M31-V1. 



M31 has been 


(Maenier et al. 


1997 


2004 


). Surveys such 



Cepheids and other variables in M31, and consequently of the underlying stellar populations 
from which they arise. What is not yet available for these variables is long temporal coverage 
with time-series data. The majority of follow-ups of Hubble's and Baa de fc S wop e's work 



are photometric in nature , focused on calibration of the Leavitt Law Jwelch et al 



1986 



Freedman fc Madorelll990l ) or on the general stellar population of M31 (iMassey et al.l 120061 ) 
rather than on the individual stars as astrophysical objects for study. This is understand- 
able given the critical importance of the Leavitt Law in modern cosmology and the obvious 
difficulty of observing stars at the distance of M31, but it gives an incomplete picture of 
the potential value of M31 Cepheids to stellar astrophysics. Measurements of period change 
in Cepheids have the potential to reveal more about the underlying stellar population and 
evolution of M31-V1 and other M31 Cepheids (ITurner. Abdel-Latif. fc Berdnikovl 120061 ) . and 
where sufficient data exist to study individual Cepheids in detail, they should be utilized. 

In 2010, observers with the American Association of Variable Star Observers (AAVSO) 
obtained just over 200 days of observations of M31-V1. These observations were in support 
of a public outreach project by the Hubble Space Telescope Hubble Heritage project, but the 
resulting data set is also a valuable new set of time series for M31-V1 in the current epoch, 
the first published time-series for this star in nearly sixty years. 

In this paper, we present our analysis of these new time series observations of M31-V1, 
and place them in context of prior observations by Hubble and by Baade. In Section 2 
we describe the observers and observations, and in Section 3 we present the results of our 
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analysis of these data. In Section 4 we present the digitized light curves of the Hubble and 
Baade & Swope data sets, and our attempt at phasing modern data with these older data 
sets. In Section 5 we conclude with a discussion of our results and make a case for future 
observations. 



Observations 



In 2010 the American Association of Variable Star Observers (AAVSO) received a re- 
quest for observations of this Cepheid by the Hubble Heritage project at STScI. Observations 
of M31-V1 were tentatively planned with the Hubbl e Space Telescope to obtain time-series 
images intended for public outreach (iNoll et al.ll201ll ). but a modern ephemeris was unavail- 
able. Such observations would be non-trivial to plan at a ground-based observatory; an 
adequate ephemeris would require at least two cycles be observed with good signal to noise, 
and a sufficient number of observations per cycle to adequately detail the light curve. For 
this star, that would require at least 60 days of coverage, with observations every 2-3 nights 
at least. For this reason, the AAVSO observer community was asked to provide coverage of 
M31-V1 for several months prior to the planned HST visit in December 2010 and January 
2011. 



M31-V1 has a magnitude around V ~ 19.4, (V — R) — 0.6 (IMassey et al.ll2006l ). which 
makes it a challenging but achievable target for amateur observers with larger telescopes and 
modern CCD cameras. The AAVSO established a set of comparison stars specifically for 
this work that lie within 15 arcminutes of M31-V1, using data from the USNO 1.0-meter and 
Sonoita Research Observatory 0.5-m telescopes. After comparison stars were established and 
made available th rough our onli ne chart plotter, the AAVSO issued a call for observations 
on July 16, 2010 (IWaagenll2010l ). Eleven different observers attempted observations of this 
target during the following six month period, and total of 214 positive observations of M31- 
VI were made between JD 2455365.9 and JD 2455565.6 (2010 June 18 to 2011 January 4). 
The resulting light curve of all AAVSO data is shown in Figured) An approximate ephemeris 
was calculated using these data, and HS T successfully im aged this region of M31 over several 
weeks, obtaining their desired data set ( iSoderbloml l201lh . 



The eleven observers who contributed data to the HST campaign are given in Table 
[TJ The majority of the observations (152 of 214 total) were obtained using Cousins R (Rc) 
filters, and our analysis is based upon these i?c~band observations. Each observer has a 
unique telescope and camera system; in general the observers were able to reach a signal 
to noise of at least 3, but a few observers were able to obtain very good signal to noise, 
with photometric accuracy much better than 0.1 magnitudes per observation. All of the 



-4- 



observations used in this paper are publicly available in the AAVSO International Database, 



and may be downloaded from the AAVSO's website (http://www.aavso.org/data-download) 
using the name "M31_V1". 

The observations are not on a common standard system, but for the purpose of time 
series analysis, this is not necessary. Observers are using common comparison stars and 
similar (if not identical) filters, which is sufficient. For our analysis we simply require that the 
amplitude each observer observes is the same; all zero-point differences were removed prior 
to time-series analysis via an iterative procedure described below. To find the relative zero 
points of each observer, we first phased the data using a rough period of 31.41 days, obtained 
by Fourier transforming the raw data from a single observer (GFB, 65 Rc observations). The 
entire data set are then folded on that period, yielding a phased light curve of all observers. 
We then divided the light curves into equal bins of 1/20 of a period, and performed an 
iterative adjustment of each observer's magnitude offset to minimize the sum of the variances 
of the 20 bins: 



EE 

i=l j=l 



N, 



o 



where iVj, is the number of bins, N Q the number of observations per bin, nij^ the magnitude 
of the j-th observation by observer k in each bin, and 8k the magnitude offset of the k-ih 
observer. Each observer's offset 8%. was iteratively adjusted up or down by a maximum of 
0.01 magnitudes per step, and the process was repeated until the sum of the variances of 
the bins reached a minimum. These offsets were then applied to each observer's data, with 
observer HQA used as the reference magnitude; the offsets are given in Table [2j and the raw 
and offset-fit phased light curves are shown in Figure [2j 



3. Analysis 



We performed a Fourier analysis on the resulting combined light curve to obtain the 
period, and phased the data to determine the time of maximum for the ephemeris. For 
the Fourier analysis, we use an iterat i ve cle aning Fourier analysis code based upon the 
algorithm of iRoberts. Lehar. fc Dreherl (119871 ) designed for irregularly-spaced observations, 
which is available from the AAVSO website. Fourier analysis of the light curve yields a peak 
at a frequency of 3.185 x 10 _2 g? — 1 (or P = 31.397c?) with an .R^-band amplitude of 0.35 
magnitudes. The precision of the frequency measurement is limited by the uncertainties 
in the observations and short temporal span of the data. The minimum 1 — a error on the 
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frequency can be estimated using the relation for a(f) given by lLenz Sz Bregerl (120051 ). which 
depends directly on the ratio of noise to amplitude, inversely on the span of the data and 
the square root of the number of data points. We assumed an average error per observation 
of 0.1 magnitudes, which yields a frequency uncertainty measured from a single season of 
data (spanning 199.8 days) of about 1 x 10 _4 c? _1 , or P = 31.4 ± O.ld. This value matches 
what was derived from earlier data within uncertainties. This is the minimum error that we 
expect on the period, and in reality it could be significantly higher if we have underestimated 
the average uncertainties of the magnitude estimates, a. 

Phasing the data using this period yields a time-of- maximum of JD 2455430.5 ±0.5 with 
the resulting linear ephemeris being 



JD max = 2455430.5(5) + 31.4(1)£. (2) 

We note that our method of establishing individual offsets may have influenced the 
determination of the period slightly, because they were calculated assuming a phasing period 
for the light curve of 31.4 days. However, we believe that this effect is minimal because the 
majority of observers obtained data at multiple phases of the light curve. If the period were 
greatly in error, it is likely that a single constant offset would not decrease the scatter across 
all phases of the light curve as we have shown. Since the bins are of order 1/20 of the period, 
we believe that the maximum error in period would be about one bin width, 1.5 days. Given 
that the period was first estimated using a single observer's data, we believe that the error 
in period is likely well below this, and is probably closer to the formal error given above. 



3.1. Color curves 

Both Rc and V band data were submitted by four different observers, three of whom 
observed this star in both filters on at least one night. We attempted to create color curves 
using data from these three observers by phasing individual observers' nightly measures of 
(V — R) with the ephemeris given above. The resulting phased color curve is shown in Figure 
[3j The phase curve is flat within the errors on individual points, but the average value of 
(V — R) is around +0.4, which is som ewhat bluer than is typical for Milky Way long period 



Cepheids (IBerdnikov fc Turner! Il995l ). However, we caution that while observers are using 
well-calibrated comparison stars for photometry, the data for these observers was not fully 
calibrated and transformed to a standard system. The fact that there are zero-point Rc- 
band differences between observers suggests that there will be V-band differences as well. 
Unfortunately there are no calibrated reference observations in V band as there are in Rc, 
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which means that our calculated average of (V — R) = 0.4 is not reliable. However, we 
can say that the variation in color over the cycle is smaller than we can measure given the 
photometric errors. 



Historic Observations: Hubble & Baade 



Both iHubbld (119291 ) and iBaade fc Swopd (119651 ) published all of their photometry in 
an easily extractable format, and we keypunched the observations from these two papers 
into electronic, machine-readble tables. These observations have since been added to the 
AAVSO International Database and are freely available via the AAVSO website. As part 
of this project, we wanted to compare how the old data compare to the new by performing 
the same analyses on the archival data a s we did on our own photom etry. In this section, 
we present analyses of the IHubbld (1l929f ) and IBaade fc Swopd (119651 ) data using the same 
procedure as in Section 3. 



4.1. Hubble photometry 



Hubble's observations were described in IHubbld (119291 ). and will be briefly reviewed 



here. The measurements were made from plates taken with the Mount Wilson 60- and 100- 
inch telescopes by ten different observers between JD 2418562.7 (1909 September 13) and 
2425149.6 (1927 September 26). We digitized 130 positive observations of M31-V1 by Hubble; 
we included all observations regardless of Hubble's quality flag ("good", "fair", or "poor") 
but we did not digitize observations where the star was below the faintest camparison, or 
where it was not visible. The resulting light curve is shown in Figure HJ 

Fourier analysis of Hubble's data yields a period of P — 31.394(i, only slightly longer 
than Hubble's own measurement of P = 31.390<i. The 1 — a error in frequency is 2.6 x 10~ 6 , 
corresponding to an error in period of ±0.003(i; our measured period and that derived by 
Hubble are identical within the uncertainties. Hubble's data folded on our period of 31.394d 
is shown in Figure El The full range amplitude of these data appears to be about the same as 
for the modern Rq data, around 1.5 magnitudes in the passband of the plates (assumed to 
be blue). It is possible that Hubble's comparison stars weren't well-calibrated, particularly 
at the faint end, and it is therefore dangerous to make assumptions about the variability by 
comparing the photographic light curve to the modern Rc curve. 
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4.2. Baade photometry 



This star was observed by Walter Baade with the Palomar 200-inch in 1950 and 1951, al- 
though the data were not published un til well after his death when Henrietta Swope brought 
his work to light ( IBaade fc Swopell965h. B aade a nd Swope's magnitudes are likely to be more 
accurate than Hubble's; IBaade fc Swopd (119651 ) used photoelectric calibration s tar m agni- 
tudes (in SA 68) from W. Baum at Palomar and lStebbins. Whitford. fc Johnson! (119501 ) . and 
Baade had much stricter control over the consistency of plates used for imaging the variable 
and calibration fields. The lightcurve of Baade & Swope's digitized data is shown in Figure 
El 

Fourier analysis of these data yields a period of 31.226d which is at first glance much 
shorter than the Hubble data. The Baade & Swope data span only two seasons (508.9 days), 
and are smaller in number than Hubble's data but the difference in period is still substantial. 
If we assume the uncertainties in magnitude are around 0.05 mag, the uncertainty in period is 
on the order of ±0.03a? which means the periods are still different by several sigma. The likely 
explanation for this discrepancy can be seen in the window function and Fourier transform 
of the Baade data, which clearly show very strong annual and monthly sidelobes (Figure 
[7]). Based upon the date s of observations, he was o bserv ing; exclusively during dark time 
around new moon. The iRoberts. Lehar. fc Dreherl (ll987| ) algorithm used for our Fourier 
calculations generally works very well in most cases, but can still produce spurious results 
when sidelobes strongly affect the transform. The problem may also be worsened given that 
the period of M31-V1 is on the order of a month which can be seen in the similarity of the 
the star's observed phase from month to month. 



Baade fc Swopd (119651 ) do not state how their period was calculated, but they appear 



to have used the Hubble period as a starting point and simply adjusted the period for a best 
fit. They give a period 31.384d for M31-V1, stating that it gives a good fit of the earlier 
Hubble data as well their own. We show a phase curve of the Baade & Swope data phased 
with our (likely incorrect period) and that obtained by S wope in Figure El Both periods 
produce acceptable phase curves, althou gh the period from IBaade fc Swopd ( 119651 ) appears 
slightly better. We therefore accept the IBaade fc Swope Jl965h period and reject our own. 
We used the 31.384d period to generate a test ephemeris with which we attempted to fit all 
measured historic and modern data together, which we discuss below. 
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4.3. Ephemeris for combined data 



As a final test, we generated an ephemeris for the entire data set from 1908 to 2010 
using the 1950-1951 da ta as a central reference point. We used the period of 31.384d from 
Baade fc Swopd ( 119651 ). along with a well-measured time of maximum of JD 2433898.0 to 
generate predicted times of maximum from 1908 to 2011. In Figure M we show these data 
overlaid onto the Rc light curve. This ephemeris produces a reasonably good fit to the 
mod ern data, fitting the max ima to within about 5% in phase. Clearly an ephemeris derived 
from iBaade fc Swopd (119651 ) can fit the modern data acceptably well, suggesting that there 
has been no measurable change in period during the entire course of observations. However, 
this cannot be guaranteed to be true, primarily because the three epochs of data are very 
widely separated, and it is possible that different periods could fit the data simply by ad- 
justing both the period and the cycle number E to produce better fits. Such a test will be 
conducted in a future paper planned once another season of data is collected. 

The primary reason that variables are observed over long periods of time or in widely- 
spaced epochs is to look for changes in variations that provide astrophys ical insight into 
their evolution. Some Cepheid variabl es are known to change over time (jSzabadoa Il983l : 
Turner. Abdel-Latif. fc Berdnikovll2006l ). and these changes can yield information abo ut their 
positio n on the H-R diagram, evol utionary state, and other physical properties. The iHubble 
( 119291 ) and lBaade fc Swopd (119651 ) data sets are a possible starting point for studying period 
changes in the M31 Cepheids. The observations by AAVSO observers represent a first 
attempt at detecting such changes serendipitously. The single season of data which we have 
obtained provides a preliminary constraint on period change in M31-V1, and we expect that 
additional seasons of data during the current epoch will provide more conclusive constraints. 
In principle, a similar modern study could be attempted on any Cepheid common to the 
Hubble and Baade & Swope samples, with the only possible limitation being the available 
telescope time and the ability of telescopes of small aperture to reach these faint magnitudes. 
We suggest that any of the M31 Cepheids with maxima brighter than R c = 19.0 would be 
suitable targets for future work. 



5. Discussion 

While the HST campaign was essentially a public outreach and education initiative, the 
resulting data set from the AAVSO community represents a new and valuable set of data for 
this star. The quality of observations suggests that the amateur community can do valuable 
work on the brighter Cepheids in M31 or nearer galaxies in the local group. The 2010 data 
set will hopefully be the first of several more years of follow-up observations. 
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The 2010 data set appears to be fully consistent with past observations of M31-V1. 
The evidence presented here - while based on only a single season of data - is consistent 
with M31-V1 being physically unchanged since its discovery by Hubble 80 years ago. It is 
certainly consistent with no significant period change since 1951, when Baade last performed 
time-series photometry, however we caution that additional seasons of data will be required 
to confirm that definitively. Subsequent years will enable us to more accurately define the 
period and times of minimum, and we encourage the continued observation of M31-V1 by the 
community, both amateur and professional. We especially encourage concentrated ifc-band 
observations around the times of maximum predicted with the Baade & Swope ephemeris. 
With well-defined times of maximum we can attempt to measure (O — C) for this star, and 
should also be able to rule out the possibility that the period has changed significantly. 

The corresponding author would like to extend his gratitude as always to the observers 
who make this research possible; the AAVSO would not exist without them. We acknowl- 
edge with thanks the variable star observations from the AAVSO International Database 
contributed by observers worldwide and used in this research. Thanks also to the AAVSO 
Sequence Team for preparation of a comparison star sequence. We also thank the referee, 
Dr. John Percy, for his very helpful comments that improved the focus of the paper. This 
research has made use of NASA's Astrophysics Data System Bibliographic Services. 
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M31-V1: All AAVSO data 




520 540 
JD - 2455000 



580 



Fig. 1. — The full AAVSO light curve of M31-V1 submitted by the eleven participating 
observers. Symbols: filled diamond, B filter; filled square, V filter; open square, clear filter 
with V zero-point; filled circle, Rc filter; open circle, clear filter with Rc zerpoint. 

Table 1. Observers contributing data to the campaign. 



Obscode 


Name 


Location 


Total 


Rc 


V c 


CR 


CV B 


GFB 


Goff 


California, US 


98 


65 


29 




4 ... 


SSTB 


Smith 


California, US 


38 


38 








HQA 


Henden 


Arizona, US 


28 


28 








SRIC 


Sabo 


Montana, US 


18 


11 


7 






BHU 


Buchheim 


California, US 


13 






13 




WGR 


Walker 


Massachusetts, US 


10 


5 


5 






HBB 


Harris 


Florida, US 


3 


2 


1 






BGU 


Belcheva 


Bulgaria 


2 




1 






CTX 


Crawford 


Oregon, US 


2 


2 








CMJA 


Cook 


Ontario, Canada 


1 


1 








DKS 


Dvorak 


Florida, US 


1 






1 
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M31-Vl:P=31.40d 
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Fig. 2. — .R-band light curves of M31-V1 before (top) and after (bottom) zero-point offsets 
were applied. Symbols: filled squares, CMJA; open squares, CTX; filled circles, GFB; 
open circles, HBB; filled triangles, HQA (reference observer); open triangles, SR1C; filled 
diamonds, SSTB; open diamonds, WGR 
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Fig. 3. — The (V — R) color curve of M31-V1 by three observers: GFB, filled circles; SRIC, 
filled squares; WGR, filled diamonds. The scatter in the color measures is large, which 
precludes us from saying anything about the color variations. 
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M31-V1: 1908-1927 
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Fig. 4. — Photographic light curve of M31-V1 from plates obtained with the Mount Wilson 
60- and 100- inch telescopes. Measurements by E. Hubble (jHubbldll929l ). 
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M31-V1: 1908-1927, P=31.394d 
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Fig. 5. — Data shown in Figure 0] folded with a period of 31.394 days. We used JD = 
2422694.3 to phase the data. 



M31-V1: 1950-1951 
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Fig. 6. — Photographic light curve of M31-V1 from plates obtained with the Palomar 200- 
inch telescope. Plates by W. Baade, photometry by Baade and H. Swope. 
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Fig. 7. — Window function (top) and Fourier transform (bottom) of the iBaade fc Swope 



(119651) data set. The window function s hows clear signs of both annual and monthly sidelobes. 



The iRoberts. Lehar. fc Dreherl (119871 ) clearly did not adequately clean the effects of the 
sidelobes, as can be seen in the double-peaked structure of the Fourier transform peaks. For 
this reason, we consider our period determination of these data to be unreliable. 

Table 2. R-band offsets for the seven observers with more than one R observation. Offsets 
were computed by minimizing the variance of the phased and binned light curve. Offsets 
are in magnitudes relative to the reference observer, HQ A. 



Obscode mag. offset N(obs) 



HQA 


0.0 


28 


CTX 


-0.10 


2 


GFB 


-0.11 


65 


HBB 


+0.27 


2 


SRIC 


+0.09 


11 


SSTB 


+0.14 


38 


WGR 


+0.29 


5 
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Fig. 8. — Phased light curve of the Baad e fc Swope data using t wo different phasing periods. 
Top: our period of 31.226 days; bottom: iBaade &: Swopel (119651 ) period of 31.384 days. Both 
periods produce acceptably good phase curves, but we believe our much shorter period is 
not correct on physical grounds that the period should not have changed by that amount 
since the early 1900s. 
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Fig. 9. — Comparison of the 2010 Rg band li ght c urve with times of maximum predicted 
using an ephemeris based on lBaade fc Swopd (119651 ). The predicted times of maximum are 
represented by the vertical lines. The predicted times are very close to the observed times, 
within 1-2 days, or about 6% of the period. 



